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Abstract 

Background: The developing vertebrate hindbrain is transiently segmented into rhonnbomeres by a process 
requiring Hox activity. Hox genes control specification of rhombomere fates, as well as the stereotypic 
differentiation of rhonnbonnere-specific neuronal populations. Accordingly, gernn line disruption of the paralog 
group 1 (PGl) /-/ox genes Hoxal and Hoxbl causes defects in hindbrain segmentation and neuron formation in 
mice. However, antisense-mediated interference with zebrafish hoxbl a and hoxbl b (analogous to murine Hoxbl 
and Hoxal, respectively) produces phenotypes that are qualitatively and quantitatively distinct from those observed 
in the mouse. This suggests that PGl Hox genes may have species-specific functions, or that anti-sense mediated 
interference may not completely inactivate /-/ox function in zebrafish. 

Results: Using zinc finger and TALEN technologies, we disrupted hoxbl a and hoxbl b in the zebrafish germ line to 
establish mutant lines for each gene. We find that zebrafish hoxbl a germ line mutants have a more severe 
phenotype than reported for Hoxbl a antisense treatment. This phenotype is similar to that observed in Hoxbl 
knock out mice, suggesting that Hoxbl/hoxbla have the same function in both species. Zebrafish hoxblb germ line 
mutants also have a more severe phenotype than reported for hoxblb antisense treatment (e.g. in the effect on 
Mauthner neuron differentiation), but this phenotype differs from that observed in Hoxal knock out mice (e.g. in 
the specification of rhombomere 5 (r5) and r6), suggesting that Hoxal /hoxblb have species-specific activities. We 
also demonstrate that Hoxblb regulates nucleosome organization at the hoxbl a promoter and that retinoic acid 
acts independently of hoxblb to activate hoxbla expression. 

Conclusions: We generated several novel germ line mutants for zebrafish hoxbla and hoxblb. Our analyses 
indicate that Hoxbl and hoxbla have comparable functions in zebrafish and mouse, suggesting a conserved 
function for these genes. In contrast, while Hoxal and hoxblb share functions in the formation of r3 and r4, they 
differ with regards to r5 and r6, where Hoxal appears to control formation of r5, but not r6, in the mouse, whereas 
hoxblb regulates formation of r6, but not r5, in zebrafish. Lastly, our data reveal independent regulation of hoxbla 
expression by retinoic acid and Hoxblb in zebrafish. 

Keywords: Zinc finger nuclease, TALEN nuclease, Retinoic acid signaling, Hindbrain, Nucleosome positioning. Gene 
expression 
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Background 

Hox genes encode a conserved family of homeodomain- 
containing transcription factors essential for metazoan 
development [1-4]. As a result of duplication events, ver- 
tebrate genomes contain four clusters of Hox genes, with 
the exception of teleost fish that have undergone an add- 
itional genome duplication - for instance, the zebrafish 
genome contains seven Hox clusters [4]. In most cases, 
genes that occupy the same position in different clusters 
(known as paralogs) have similar expression patterns 
and functions, leading to redundancy of Hox function. 
During early development, Hox genes specify tissue 
identities along the anterior-posterior (AP) axis of the 
animal. The linear arrangement of Hox genes in the gen- 
omic clusters coincides with the timing and position of 
their expression along the AP axis, a characteristic 
termed colinearity [3,5,6]. The retinoic acid (RA) signal- 
ing pathway activates early Hox gene expression and is 
important in colinear regulation [7,8]. RA binds a het- 
erodimeric complex of RA receptors (RARs) and retinoic 
X receptors (RXRs) that target cis-regulatory sites 
known as RA response elements (RAREs) in the Hox 
clusters [9-11]. RA promotes decondensation of Hox 
clusters from compact chromosomal chromatin in cells 
and embryos [12-14] and this process correlates with the 
progressive activation of Hox transcription along the 
genomic cluster. Once transcribed, Hox genes also regu- 
late the expression of other Hox genes in auto- and 
cross -regulatory loops. The highly conserved process of 
Hox gene activation and regulation leads to an overlap- 
ping series of Hox expression domains along the AP 
axis, sometimes referred to as the "Hox code" [15]. 

During early embryogenesis, the presumptive verte- 
brate hindbrain is transiently divided into seven to eight 
segments (rhombomeres) and Hox genes play a key role 
in formation of the more posterior rhombomeres [16]. 
Each rhombomere gives rise to unique cell populations 
from which segment-specific motor neurons and reticu- 
lospinal neurons differentiate. For motor neurons, this 
includes the trigeminal neurons in rhombomere 2 (r2) 
and r3, the facial motor neurons (FMNs) in r4, the abdu- 
cens neurons in r5 and r6, and the vagal neurons in the 
caudal hindbrain. These neuronal pools form the motor 
neuron nuclei of the V'^, Vl'^, VIl'^ and X'^ cranial 
nerves that innervate the face, head, and neck of the ani- 
mal. Reticulospinal neurons also form in a rhombomere 
determinate manner and regulate breathing and circula- 
tion, as well as the coordination of locomotor signals 
between the spinal cord and the brain. For instance, the 
Mauthner neurons (MNs) form in r4 and extend axons 
contralaterally and posteriorly and function to control 
the escape response. Segmentation of the hindbrain 
starts with the formation of r4 followed by rl/r2, r3, r7, 
and r5/r6 [17]. Accordingly, the first Hox genes transcribed 



in the mouse, Hoxal and Hoxbl, are expressed in r4 
[5,18,19]. A series of loss-of-function studies have deter- 
mined that Hoxal and its downstream target Hoxbl 
have separate functions. In particular, Hoxal mutant 
mouse embryos have segmentation defects while Hoxbl 
mutants appear to have neuronal defects related to r4 
specification. The segmentation defects observed in 
Hoxal mutants include an enlarged r3, a reduced r4, 
and a reduced or completely lost r5 [20-24]. Similar seg- 
mentation defects are also found in mice with mutations 
made to the retinoic acid response element (RARE) 
found in the downstream enhancer of Hoxal [25] . These 
segmentation defects are specific to the function of 
Hoxal, as mouse Hoxbl mutant embryos show no defects 
in hindbrain segmentation [26,27]. While hindbrain seg- 
ments form normally in Hoxbl mutants, r4 derived FMNs 
fail to migrate into r5 [23,26,28] and instead migrate away 
from the midline of r4, assuming lateral positions similar 
to trigeminal neurons in r3. FMN neurons in Hoxbl mu- 
tants extend axons out through r4 into the second 
pharyngeal arch similar to the projection of wild type 
FMN axons. Loss-of-function studies for zebrafish hoxblb 
and hoxbl a (functionally analogous to murine Hoxal and 
Hoxbl, respectively) have been accomplished using anti- 
sense morpholino oligos (MOs) to block translation of 
hoxblb and hoxbl a. Embryos injected with hoxbl bMO 
have hindbrain segmentation defects with an expanded r3 
and a reduced r4, r5 and r6, while hindbrain segmentation 
is unaffected in hoxbl aMO-injected embryos [29]. Fur- 
thermore, hoxbl aMO'injected embryos possess FMNs 
remaining in r4 that resemble the stalled FMNs observed 
in mouse Hoxbl mutants. While these data indicate that 
the zebrafish hoxblb and hoxbla genes have roles similar 
to those of mouse Hoxal and Hoxbl, there are also differ- 
ences between the defects observed in the mouse mutants 
versus MO-injected zebrafish. First, segmentation defects 
appear more severe in Hoxar""' and Hoxar^~ ;Hoxbr^~ 
mutants than in hoxbl bMO and hoxbl bMO; hoxbl aMO- 
injected zebrafish embryos. Specifically, the r5 domain in 
Hoxar^~ mouse embryos is lost, while r5 is merely reduced 
in zebrafish hoxbl bMO-injected embryos. Hoxal~^~; 
Hoxbr^~ mice also appear to have a stronger segmenta- 
tion defect, with loss of both r4 and r5, while hoxbl bMO; 
hoxbl aMO'injected zebrafish embryos show only a 50% 
reduction in the size of r4. Second, unlike the HoxaF^' 
mice, hoxbl bMO'injected embryos have a reduced r6. 
Third, the hoxbl aMO shows no effect on the reticulosp- 
inal neurons in r4 although Hoxbl mutant mice have a 
miss-specified r4. Indeed, even when hoxbl bMO and 
hoxblaMO are co-injected, only incomplete formation of 
Mauthner neurons is observed in r4. These phenotypic 
differences suggest that either these genes function dif- 
ferently in mouse and zebrafish, or the MO phenotypes 
do not recapitulate true loss-of-function of hoxblb and 
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hoxbla. Here we report the generation of hoxbla and 
hoxblb germ line mutations using zinc finger and TALE 
nucleases in zebrafish. We find that these germ line mu- 
tants have more severe phenotypes than those reported 
from antisense MO injections, suggesting that the MOs 
do not completely block Hox function. Our results dem- 
onstrate that hoxbla is required for specification of r4 
and for the formation of r4-specific neurons. This is 
similar to the reported mouse phenotype, suggesting 
that hoxbla and Hoxbl share conserved functions in r4 
formation. Similarly, our results reveal that zebrafish 
hoxblb and mouse Hoxal share functions in the forma- 
tion of r3 and r4, but their roles also differ, such that 
Hoxal is required in mouse r5 while hoxblb is required 
in zebrafish r6. Lastly, we demonstrate that hoxblb and 
retinoic acid act separately to activate hoxbla expres- 
sion, with Hoxblb modulating nucleosome organization 
at the hoxbla promoter. 

Results 

Generation of hoxbla and hoxblb germ line mutants 

To investigate the roles of hoxbla and hoxblb in zebra- 
fish hindbrain development, we generated hoxbla and 
hoxblb loss of function mutants using zinc finger nucle- 
ases (ZFNs) and TALE nucleases (TALENs). ZFNs and 
TALENs consist of the Fokl endonuclease tethered to 
sequence-specific DNA-binding domains (zinc finger or 
TALE) that target the nuclease to a desired genomic lo- 
cation [30,31]. The use of obligate heterodimeric Fold 
nucleases increases target specificity by requiring that 
two complementary Fold nuclease domains bind to 
adjacent genomic sequences for activity [32,33], thereby 



reducing, but not eliminating, the likelihood of off-target 
effects. Once activated through dimerization, the Fokl 
nuclease introduces a double strand DNA break that is 
repaired through the non-homologous end joining (NHEJ) 
repair pathway. NHEJ is error prone and will introduce 
mutations at a low rate. XXHiile many of the resulting mu- 
tations do not affect protein function, we were particularly 
interested in identifying the small number of mutations 
that lead to shifts in the reading frame and introduce stop 
codons. 

We initially employed ZFNs to target both hoxbla and 
hoxblb based on several criteria. First, nucleases were 
targeted to a site in the first exon of each gene in order 
to increase the likelihood that a frame shift would ter- 
minate translation upstream of known functional do- 
mains, particularly the homeodomain. Second, the spacing 
between the target sequences for each ZFN pair was set to 
either 5 bp or 6 bp based on previous reports indicating 
that these represent optimal spacing [34]. Third, we tar- 
geted regions containing a restriction site that could be 
used to screen for mutations. Based on these criteria, we 
designed several ZFNs to each gene using three separate 
ZFN "builds" (Table 1). The first and second builds were 
based on a modular library of single zinc finger proteins 
[35] . For build 1, we generated one ZFN targeting hoxbla 
(Zbla-1) and one targeting hoxblb (Zblb-1) where each 
ZFN contained three zinc fingers (ZFs) assembled from 
the modular library. Build 2 (Zbla-2 and Zblb-2) 
employed the same modular library and targeted the 
same genomic sites as build 1, but each fusion protein 
contained four ZFs to increase specificity for the target 
sequence while decreasing the instances of off-target 



Table 1 Characteristics of TALE and zinc finger nucleases targeting hoxbla and hoxblb 



TALEN/ZFN 


Target 
coordinate^ 


Upstream TALEN/ZFN 
Target sequence'^ 


Downstream TALEN/ZFN 
Target sequence'^ 


Gap" 


No. of units^ 


Embryos 
injected^ 


Activity^ 


Diagnostic 
enzyme^ 


Zblb-1 


Chrl 2:2871 2770 


GGTGGAAGG 


GTGGACATG 


5bp 


3 


274 


NO 


Bsll 


Zblb-2 


Chrl 2:2871 2770 


GGTGGAAGGGCT 


GTGGACATGGGT 


5bp 


4 


424 


NO 


Bsll 


Zblb-3 


Chrl 2:2871 2770 


GTGGAAGGGCTG 


GTGGACATGGGT 


6bp 


4 


109 


YES 


Bsll 


Zblb-4 


Chrl 2:2871 2773 


GAGGTGGAAGGG 


GACATGGGTAAA 


6bp 


4 


149 


YES 


Bsll 


Zbla-1 


Chr3:24060660 


GCTGATAAG 


GATGCGAAG 


6bp 


3 


2916 


NO 


Rati 


Zbla-2 


Chr3:24060660 


GCTGATAAGATG 


GATGCGAAGGCC 


6bp 


4 


541 


NO 


Rati 


Zbla-3 


Chr3:24060602 


GCCATAGTGTGG 


GCCGGTGCGTAC 


6bp 


4 


72 


NO 


Bsll 


Zbla-4 


Chr3:24060535 


AGGGTOATAAA 


GGATGGGATGTA 


5 P 


4 


56 


NO 


Bsll 


Tbla-1 


Chr3:24060209 


TCCAGAATGAACTC 


TCCCACGGTOCAAAT 


16bp 


15/16 


50 


NO 


Rsal 


Tbla-2 


Chr3:24060227 


CTOGAGTACACAAT 


TGGGCGAGTAGGCG^ 


16bp 


16/16 


50 


YES 


Btgl 


Tbla-3 


Chr3:24060213 


CCAGAATGAACTCmC 


TCGTCCCACGGTOC 


16bp 


18/15 


50 


NO 


Rsal 



^Target coordinate as defined by midpoint between the upstream and downstream nuclease recognition sites. 

"^Genomic sequence targeted by the upstream and downstream TALE and ZF nucleases. 

^Distance between upstream and downstream TALE and ZF nucleases. 

^The number of TAL or ZF units used to assemble each nuclease. 

^The number of embryos injected to test activity and raise mutant families. 

^Nucleases that disrupted the diagnostic restriction site were defined as active. 

^The enzyme used to test activity of nucleases. 
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effects. For build 3, we designed two ZFNs to each gene 
(Zbla-3, Zbla-4, Zblb-3 and Zblb-4). Zblb-3 and 
Zblb-4 targeted the same genomic sequence as Zblb-1 
and Zblb-2, although Zblb-4 was offset 3 bp relative to 
Zblb-3, while Zbla-3 and Zbla-4 targeted sites 60 bp 
and 125 bp, respectively, upstream of the site targeted 
by Zbla-1 and Zbla-2. The ZFNs designed in build 3 also 
used four ZFs each, but were assembled from an updated 
version of the ZF library that includes two-finger modules 
[36]. These two-finger modules were optimized for their 
ability to bind DNA efficiently in tandem to minimize the 
impact of context-dependent effects on recognition when 
fingers are assembled into arrays. 

In vitro transcribed mRNA encoding each ZFN pair 
was injected into early one-cell stage embryos and gen- 
omic DNA was prepared from pools of whole embryos 
collected 24 hours post fertilization (hpf; Figure lA). 
ZFN activity was measured by amplifying the targeted 
region, followed by digestion to estimate the fraction of 
genome with a disruption of the diagnostic restriction 
site. Notably, the ZFNs are likely to act after the first 
several cell divisions (due to the rapid cell cycle of zebra- 
fish embryos, as well as due to the need for the ZFN 
mRNA to be translated) and it is therefore expected that 
mutations will be induced in only a subset of cells - ren- 
dering the embryos mosaic. Embryos injected with 
Zblb-3 and Zblb-4, but not Zblb-2 or Zblb-1 (Figure IB) 
revealed a partial loss of the diagnostic restriction site, 
suggesting that these ZFNs are active. Based on the inten- 
sity of the uncut band in the diagnostic digest, the Zblb-3 
ZFN may be more active than Zblb-4. Since ZFNs from 
all three builds failed to induce mutations at the hoxbla 
locus, we turned to TALENs as an alternative method to 
disrupt the hoxbla gene. To increase the likelihood of 
success, we generated three different hoxbla TALENs that 
differ slightly in the length of their target sequences 
(Table 1) using Golden Gate TALEN assembly [31]. As 
with ZFNs, TALENs were designed to target regions in 
the first exon of hoxbla that include a diagnostic restric- 
tion site. Notably, the TALENs were directed to sites 
450 bp (Tbla-1) or 430 bp (Tbla-2 and Tbla-3) upstream 
from the region targeted by hoxbla ZFNs. Using the same 
mRNA microinjection strategy as for the ZFNs, we found 
that TALEN Tbla-2 introduced mutations - as evidenced 
by loss of the diagnostic restriction site - but that Tbla-1 
and Tbla-3 did not (Figure IB). 

Having identified functional ZFNs and TALENs, we 
raised injected embryos to establish an adult FO founder 
population. As noted, fish in this FO population will be 
mosaic and each individual fish may carry more than 
one mutant allele for the same gene (Figure lA). Geno- 
typing of pooled offspring identified 20 hoxbla (out of 
24 tested) and 15 hoxblb (out of 35 tested) FO founders 
that transmit mutations via their germ lines, suggesting 



mutagenesis rates of 83% and 43% for hoxbla and 
hoxblb, respectively. Sequencing of mutant alleles from 
each FO founder revealed that two of the 20 hoxbla (A2 
and A20) and three of the 15 hoxblb (B2, Bll and B15) 
FO founders carry mutations that introduce frame shifts 
(Table 2), while the remaining FO founders transmitted 
mutations that disrupt the diagnostic restriction site, but 
that do not create frame shifts (two representative exam- 
ples are shown in Additional file 1: Table SI). The five 
FO fish that carry mutations causing frame shifts were 
outcrossed to wild type fish and the resulting offspring 
raised to generate the Fl generation. Genotyping of Fl 
fish allowed us to determine the transmission rate of 
mutations from mosaic FO fish. We find that the hoxblb 
FO founders transmit their mutations at a frequency of 
--40% (45% for B2, 43% for Bll and 41% for B15), while 
the hoxbla FO fish transmit their mutations at -10% 
(14% for A2 and 9% for A20). Notably, we also find that 
the three hoxblb founders each transmitted only one 
mutant allele, while the two hoxbla founders transmit- 
ted multiple mutant alleles each (four alleles from A2 
and two alleles from A20). Thus, we have generated six 
hoxbla {uml89, uml90, uml91, uml92, uml93, uml94) 
and three hoxblb (uml9S, uml96, uml97) mutant lines 
(Table 2) that were used for further phenotypic analysis. 

Closer analysis of the mutant sequences revealed that 
Zblb-3 and Tbla-2 generated both insertions and dele- 
tions (Figure IC). In particular, the Zblb-3 ZFN intro- 
duced deletions ranging from 1 bp to 8 bp, as well as a 
19 bp insertion in the hoxblb gene, while the Tbla-2 
TALEN introduced deletions ranging from 1 bp to 8 bp, 
as well as a 17 bp insertion, in the hoxbla gene. We note 
that large insertions and deletions that interfere with the 
PGR reaction (e.g. by deleting a primer site) would not 
be detected by our experiments, suggesting that the sizes 
observed here may be somewhat biased to smaller dele- 
tions and insertions. Conceptual translation of each mu- 
tant allele confirmed a shift in the reading frame 
(Figure 1 and Additional file 2: Figure SI). As a result, 
hoxbla mutant alleles go out of frame starting with resi- 
due 15 {uml93)y residue 16 {uml90, uml91 and uml92) 
or residue 17 (uml89 and uml94) and hoxblb mutants 
starting with residue 71 {uml97), residue 74 {uml9S) or 
residue 75 {uml96). While the out of frame sequences 
code for varying numbers of missense residues, all ter- 
minate in a premature stop codon and none of the mu- 
tant alleles is predicted to encode a homeodomain. 

Lastly, we raised embryos from hoxbla^^^^^^^ x 
hoxbla^^"""^'^^ and hoxblb^ ^"""^'^^ x hoxblb^ ^"""^'^^ 
crosses to adulthood and genotyped them. We find that 
hoxbW'^^^^^'''^^^^ animals do not survive to adulthood 
(0/23 genotyped adults were hoxbla'''^'^'^'''^'^^), while 
hoxblb'''^'^^^'''^'^^ embryos do - although their viability 
may be somewhat reduced (13/84 genotyped adults were 
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1 . Inject wildtype embryos at one cell stage 



hoxbia 



I 

- Screen embryo pools by restriction 
digest to confirm ZFN/TALEN activity 

2. Raise mosaic founders (F0) 



25bp 



A 



- Screen embryo pools by restriction digest 
to confirm germline transmission of lesion 
- Sequence lesion to identify frameshift 

3. Raise carriers (F1) 



Btgl site 



bl a"^ TTTGTAACCGTGGGACG AACGCC 

^^QumisP JTJGTAAC GTGGGACGAACGCC -Ibp 

^lQumi90 TTTGTACACAATTTGTACAATTT +17bp 

^^^amiPi jjTGTAC GGACGAACGCC -5bp 

^ -J ^umi92 J J jGTACTCCATTTGTA C -5bp 

t)ia^n,i93 jj TGGGACGAACGCC -8bp 

bla"'"^^^ TTTGTAATTT GGACGAACGCC -2bp 




- Screen for lesion to identify carriers 
- Sequence lesion of carriers to confirm 
lesions that create frameshift 

4. F1 carriers crossed for function analysis 

X 



bla^ 

|3l3uml89 I 
|33_3uml91 

203bp 



!l7 aa 

!l6 aa 

■ 

!l6 aa 

■ 

'l6 a a 
■ 

is aa 
■ 

!l7 aa 



- 316 aa 

- 55 aa 

- 59 aa 

- 42 aa 

- 42 aa 

- 41 aa 

- 43 aa 



hoxbib 



F1 



* 



F1 




+/+ +/- 




3 



Bsll^site 

blb"^' ATCGCAGCCCTTCCACCTCAACGTGGACA 
^^^.mi95 ATCGCAGCCCTTCC GTGGACA -8bp 

^Ib"'"^^^ ATCGCAGCCCTTCCACATTCC GTGGACA -Ibp 

blb"""^' ATCGTGTTTCAATCGTGAAACACAAATTCACAAATTCAATCGTG +19bp 



blb^^ 








aa 








- 77 


aa 




V3 


aa 












- 120 


aa 




V4 


aa 












- 86 


aa 




Ve 


aa 







F2 

B > > 

^ ^ ^ <<y 

Figure 1 Generation of hoxbia and hoxbib germ line mutants using ZFN and TALEN technologies. A. Diagram outlining experimental 
strategy. Active ZFNs and TALENs were identified by their ability to disrupt the sequence of a diagnostic restriction site in genomic DNA. Embryos 
injected with active ZFNs and TALENs were then raised and screened for founders that transmit frameshift mutations via their germ line. B. 
Identification of active nucleases. Genomic DNA was prepared from pools of injected embryos and digested with a diagnostic restriction enzyme. 
The Zblb-3, Zblb-4 and Tbla-2 injected pools contain undigested material (arrows), indicating that the diagnostic restriction site has been 
disrupted. C. Characterization of germ line transmitted hoxbia and hoxbib mutations. Six mutant hoxbia (top) and three mutant hoxbib (bottom) 
alleles were identified to cause frameshift mutations. Nucleotide and peptide sequences are indicated for each mutant allele. Gaps in the 
nucleotide sequence indicate deletions, while red nucleotides indicate insertions. Numbers to the right of each nucleotide sequence indicate the 
net size of insertions/deletions. For the peptide sequence, gray boxes indicate residues read out of frame prior to encountering a premature stop 
codon. Amino acid numbers below each peptide sequence indicate the residue affected by the frameshift mutation, while numbers to the right 
indicate the length of the mutant peptide. HD = homeodomain. Black wedge indicates site of single intron in each sequence. 
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Table 2 Characteristics of hoxbia and hoxbib germ line mutant alleles 


Founder 


Transmission frequency 


Sequence^ 


Type of mutation'^ 


ID^ 


hoxbia 




mGTAACCGTGGGACGAACGCCTACT 






A2 


14% 


mGTAAC GTGGGACGAACGCCTACT 


Deletion (-1 bp) 


uml89 






mGTAC/\C/\/\777G7/\C/\/\777GGACGAACGCCTACT 


Insertion/deletion (+17 bp) 


urn 190 






JTTGJACGGG GGACGAACGCCTACT 


Insertion/deletion (-5 bp) 


uml91 






TTTGJACTCCATTTGTA CTACT 


Insertion/deletion (-5 bp) 


urn 192 


A20 


9% 


TT TGGGACGAACGCCTACT 


Deletion (-8 bp) 


uml93 






mGTAATTTC GGACGAACGCCTACT 


Insertion/deletion (-2 bp) 


urn 194 


hoxbib 




ATCGCAGCCCTOCACCTCAACGTGGACATGGG 






B2 


45% 


ATCGCAGCCCTOCACCTCAACGTGGACATGGG 


Deletion (-8 bp) 


uml95 


B11 


43% 


ATCGCAGCCCTOCAC/\7"7"CC GTGGACATGGG 


Insertion/deletion (-1 bp) 


urn 196 


B15 


41% 


7G7"77"C/\/\7CG7"G/\/\/\C/\C/\/\/\7"7"C/\C/\/\/\77"C/\/\7"CGTGGACATGG 


Insertion/deletion (+19 bp) 


urn 197 



^Sequence of mutant alleles. Deletions are shown as gaps and insertions are italicized. Nucleotides in bold indicate 
wild type sequence. 

"^Indicates whether mutation results from insertion, deletion or both. Numbers in parenthesis indicate net gain/loss 
'^Identifying designation for each mutant allele. 



target site for diagnostic restriction enzyme in 
of nucleotides in mutant sequence. 



hoxblb'''^'^^^'''^'^^). We note that this is in contrast to 
the situation in the mouse, where both Hoxal'^' and 
Hoxbr^~ pups die shortly after birth [20,26]. 

hoxbib is required for zebrafish hindbrain segmentation 

Formation of the vertebrate hindbrain requires segmen- 
tation of the neural tube into rhombomere domains, 
as well as the specification of distinct cell fates and the 
differentiation of characteristic types of neurons in each 
rhombomere. Paralog group 1 (PGl) hox genes, such 
as hoxbia and hoxbib, are among the earliest genes 
expressed in the hindbrain primordium and hox function 
has been implicated in multiple aspects of hindbrain de- 
velopment (reviewed in [16]). We therefore made use of 
the hoxbia and hoxbib mutant lines to examine the role 
of these genes in development of the zebrafish hindbrain. 

We first examined the expression of several rhombomere- 
restricted genes - pax2 (expressed at the midbrain- 
hindbrain boundary; MHB), krox20 (expressed in r3 and 
r5), hoxbia (expressed in r4), hoxbSa (expressed in r5 
and r6) and hoxd4a (expressed in r7 and r8). For this 
purpose, heterozygous hoxbib^ ^^"^^^^ Fl fish were in- 
crossed and the resulting embryos were assayed at 
22hpf by in situ hybridization followed by genotyping. 

mutant 
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We find that homozygous hoxbW 
embryos express krox20 in r3 and r5, as well as hoxbia 
in r4 (Figure 2C). However, the size of r3 is increased 
and the size of r4 is decreased in hoxbib mutants rela- 
tive to wild type (or heterozygous) embryos (Figure 2A; 
Additional file 3: Table S2). To address the possibility 
that the Zblb-3 ZFN might have introduced off-target 
mutations in the hoxblb^^^^^ line that could contribute 
to this phenotype, we also examined in-crosses of the 
hoxblb^^^^^ and hoxblb^^^^^ lines, as well as pair-wise 
inter-crosses among all three lines. We find that mutant 



embryos derived from all such crosses exhibit the same 
phenotype and that this phenotype segregates with the 
hoxbib mutation (Additional file 3: Table S2 and Additional 
file 4: Figure S2), suggesting that it is due to disruption of 
the hoxbib gene. Furthermore, the PROGNOS on-line tool 
[37] revealed five exonic sites in the top fifty potential 
off-target sites for Zblb-3, but neither of these sites re- 
sides on the same chromosome as the hoxbib gene, and 
they would therefore segregate independently of the 
hoxbib mutation in our crosses. Further analysis of 
hoxblb^^^^^^^^^^^ mutant embryos revealed expression 
of pax2, hoxbSa and hoxd4a in the expected domains 
(Figure 2G, O). In addition to the enlargement of r3 and 
the reduction of r4 noted above, this analysis also re- 
vealed an apparent reduction of r6 - as evidenced by a 
smaller gap between r5 krox20 staining and r7 hoxd4a 
staining (brackets in Figure 2E, G), as well as by a re- 
duction in the size of the hoxbSa expression domain 
(brackets in Figure 2M, O) in mutant embryos relative 
to wild type embryos. 

We next quantified the changes in rhombomere size 
by direct measurements (Figure 2P). We find that r3 is 
significantly enlarged (79.4 (im in hoxbib mutants versus 
62.0 (im in wild type; p = 0.0003) and r4 significantly re- 
duced (18.5 (im in mutant versus 47.4 (im in wild type; 
p = 2.53E-12) in mutant embryos. Notably, we cannot dis- 
tinguish whether this effect is due to some cells switching 
from an r4 to an r3 fate, or if r3 cells have a growth ad- 
vantage in the absence of hoxbib function. We also find 
that r6 (26.5 (im in mutant versus 47.6 (im in wild type; 
p = 3.32E-09) and rl/r2 (129 (im in mutant versus 
139 (im in wild type; p = 0.007) are somewhat reduced, 
but r5 is unaffected, in hoxbib mutants. Accordingly, 
measuring the length of the entire hindbrain reveals it 
to be significantly shorter in hoxbib mutant embryos 
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Figure 2 Zebrafish hoxblb is required for hindbrain 
segmentation. A-0. Wild type (A, E, I, M), hoxblo~^~ {B, f, J, N), 

hoxblb~^~ [C, G, K, O) or doubly hoxblQ~^~;hoxblb~^~ embryos 
(D, H, L) were assayed by in situ hybridization for expression of 
hoxblo in r4 (blue stain in A-D), kfox20 in r3/r5 (red stain in A-D, 
1-0 and blue stain in panels E-H), pQx2 at the midbrain-hindbrain 
boundary (blue stain in E-H), fgfS in r4 (blue stain in l-L), hoxd4o in 
r7 (blue stain in E-H) and hoxbSo in r5/r6 (blue stain in M-0). 
P. Quantification of segmentation defects in hoxblb~^~ embryos. 
Rhombomere lengths were measured as indicated in the inset. 
MHB-r6 measures the full distance from the anterior limit of the 
MHB to the posterior limit of r6. p-values were computed using 
Students' t-test and error bars represent standard error. N = 10 
embryos. All embryos are flat mounted in dorsal view with anterior 
to the top. A-H and M-0 are at 22hpf, while l-L are at 14hpf 
r = rhombomere; MHB = midbrain/hindbrain boundary. 



(306 \im in mutant versus 354 \im in wild type; p = 
3.28E-05), presumably as a result of the reduced length 
of several rhombomeres. 

Analysis of hoxbla mutant embryos (derived from an 
inter-cross of hoxhld''^'^'^^ and hoxhla!''^'^^^^ carriers) 
revealed normal expression of krox20 in r3 and r5 of the 
hindbrain (Figure 2B, F, J, N). Furthermore, the size of 
the r4 domain is normal, but expression of hoxbla is 
markedly reduced in r4 of hoxbla mutants (Figure 2B; 
Additional file 3: Table S2). Examination of embryos de- 
rived from inter-crosses of the hoxbla^^^^^, hoxbla^^^^^, 
hoxbla'''^'^^ and hoxbla"""^'^^ alleles confirmed this 
phenotype and suggest that it results from mutation of 
hoxbla rather than from off-target mutations introduced 
by the Tbla-2 TALEN (Additional file 3: Table S2 and 
Additional file 5: Figure S3). Accordingly, the PROGNOS 
tool did not reveal any exonic sites among the top fifty po- 
tential off-target sites for the Tbla-2 TALEN. Since 
hoxbla regulates its own expression [38], the loss of 
hoxbla transcript in hoxbla mutants may be due to the 
loss of Hoxbla protein, or, alternatively, to reduced stabil- 
ity of the mutant hoxbla transcript. Accordingly, expres- 
sion of fgf3, another r4-restricted gene dependent on 
hoxbla for expression, is also markedly reduced in r4 of 
hoxbla mutants at 14hpf (Figure 2J). Further analyses 
demonstrated normal expression of pax2, hoxd4a and 
hoxbSa (Figure 2F, N) and also revealed normal size of 
rhombomeres in hoxbla mutant embryos. We conclude 
that hoxblb is required for formation of appropriately 
sized rhombomere segments in the caudal hindbrain. In 
contrast, hoxbla is essential in regulating expression of 
r4-restricted genes, indicating a role for hoxbla in r4 
specification. 

hoxbla is required for formation of r4-derived neurons 

A key event in hindbrain development is the differenti- 
ation of unique complements of neurons in each rhombo- 
mere. In particular, motor neurons of the V^^ (trigeminal) 
cranial nerve differentiate in r2 and r3, motor neurons of 
the Vf^ (abducens) cranial nerve form in r5 and r6 and 
motor neurons of the X^^ (vagal) cranial nerve form in the 
caudal most region of the hindbrain. In addition, motor 
neurons of the VII^^ (facial) cranial nerve form in r4, but 
subsequently migrate to r6 and r7 in zebrafish. In order 
to determine if neuronal differentiation is affected in 
hoxbla and hoxblb mutants, we analyzed cranial motor 
neuron formation. Wild type embryos displayed the 
expected stereotypical arrangement of cranial motor 
neurons at 48hpf (Figure 3A). Notably, this includes an 
almost complete lack of facial motor neurons (FMNs) in 
r4 as a result of these neurons having migrated caudally 
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mutant 



by this stage. In contrast, all hoxbla 
embryos revealed a large number of motor neurons in 
r4 and reduced numbers in r6 and r7 (Figure 3B). We 
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Figure 3 hoxbla and hoxblb are required for neuronal differentiation in the hindbrain. 48hpf wild type (A, F), hoxbla~^~ {B, G), hoxblb 

(C D, H, I) or doubly hoxbla~^~:hoxbl b~^~ embryos (E, J) were assayed by immunostaining for the differentiation of branchiomotor neurons 

(isletl/2 staining in A-E) and Mauthner neurons (3A10 staining in F-J). Colored arrowheads indicate r4 (yellow), r5 (blue) and r6 (red). All embryos 

are flat mounted in dorsal view with anterior to the top. 
I. J 



interpret this to mean that FMNs fail to migrate out of 
r4 and that the residual motor neurons in r5/r6 primar- 
ily represent abducens neurons. FMNs also remain in r4 
of all hoxblb mutant embryos, but the phenotype is 
more variable than what is observed in hoxbla mutants. 
Approximately half of the hoxblb mutants reveal sig- 
nificant retention of FMNs in r4 (Figure 3C), while the 
other half shows more extensive FMN migration (Fig- 
ure 3D). In addition to impaired FMN migration from 
r4, the number of FMNs appears lower in hoxblb mu- 
tants and these are less well organized, with cells being 
less tightly grouped in r5-r7 and numerous cells found 
located outside the main clusters. 

Similar to the cranial motor neurons, reticulospinal 
neurons also display rhombomere-specific differentiation 
in zebrafish. Specifically, bilaterally arranged Mauthner 
neurons form in r4 and project their axons across the 
midline down into the spinal cord at 36hpf (Figure 3F). 
Notably, differentiated Mauthner neurons are absent in 
hoxbla mutant embryos (Figure 3G). As observed for 
the cranial motor neurons, hoxblb mutants show a vari- 
able phenotype such that approximately half of the em- 
bryos retain one Mauthner neuron, while the other half 
fails to form both Mauthner neurons (Figure 3H, I). We 
conclude that hoxbla function is absolutely required for 
FMN migration and Mauthner neuron formation, but 
that hoxblb is only partially required for these processes. 

hoxbla and hoxblb have separate functions in zebrafish 
hindbrain development 

Since hoxbla and hoxblb are both required for normal r4 
formation, we examined their functional relationship by 
analyzing hoxbla/hoxblb double mutant embryos (gener- 
ated from a cross between hoxbla^^"""^'^^; hoxblb^ ^"""^'^^ 



and hoxbla^''''^'^^; hoxblb^ ''''^'^^ double heterozygotes). 
Using in situ hybridization, we find that double mutant 
embryos have hindbrain segmentation defects with an ex- 
pansion of r3 and a reduction of r4 and r6 (Figure 2D, H; 
Additional file 3: Table S2). These changes in rhombo- 
mere size are indistinguishable from those observed in 
hoxblb mutant embryos (compare Figure 2G to H), fur- 
ther demonstrating that hoxbla does not play a role in 
zebrafish hindbrain segmentation. Furthermore, double 
mutant embryos lack hoxbla expression, while single 
mutants show a reduction either in the level (hoxbla 
mutants) or domain size (hoxblb mutants) of hoxbla 
expression (compare Figures 2B-D), indicating that both 
genes may play a role in hoxbla transcription. 

Next we examined neuronal differentiation in double 
mutant embryos. We find that FMNs form in r4 of 
double mutants, but do not migrate caudally (Figure 3E), 
similar to the phenotype of hoxbla mutants (Figure 3B). 
We also note that the population of FMNs is smaller in 
double mutants (Figure 3E), similar to the phenotype 
of hoxblb mutants (Figure 3C, D). Lastly, double mu- 
tants completely lack Mauthner neuron formation in r4 
(Figure 3J), identical to the hoxbla mutant phenotype. 
These findings are consistent with hoxblb being required 
for the size of r4 (and therefore the number of neu- 
rons formed in r4), while hoxbla is required for the 
migration of FMNs and the differentiation of Mauthner 
neurons. Notably, the abducens neuron population in 
r6 may be slightly smaller in double mutants, consist- 
ent with r6 being smaller in the absence of hoxblb, but 
this effect cannot be seen in hoxblb single mutants due 
to the residual migrating FMNs. We conclude that 
hoxbla and hoxblb have different functions in hind- 
brain development. 
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hoxbib and retinoic acid independently activate hoxbia 
expression 

Although the size of r4 is reduced in hoxbib mutant em- 
bryos, we note that hoxbia expression persists. The ex- 
pression of hoxbia in the absence of hoxbib is somewhat 
surprising since previous work suggested that early neural 
expression of hoxbia depends on a Hoxb lb-regulated en- 
hancer [29]. Since retinoic acid (RA) is known to activate 
hox gene transcription in many settings [13,25,28,39-43], 
we next investigated whether RA signaling activates 
hoxbia transcription in the absence of hoxbib. 

To this end, embryos were treated with lOuM diethy- 
laminobenzaldehyde (DEAB; a small molecule inhibitor 
of the RALDH enzyme involved in RA synthesis), or 
with lOOnM exogenous retinoic acid, from Ihpf to 19hpf 
(Figure 4). We find that, in wild type embryos, lOuM 
DEAB blocks krox20 expression in r5, though not in r3 
(that instead appears to be expanded), but does not 
affect hoxbia expression in r4 (Figure 4C). Treatment of 
wild type embryos with RA produces a distinct phenotype 
such that hoxbia expression is expanded, while krox20 ex- 
pression is lost in r3 and reduced in r5 (Figure 4E). These 
results indicate that RA supports krox20 expression in r5 
while it inhibits it in r3, consistent with previous reports 
that RA promotes posterior fates [44,45]. Furthermore, 
while hoxbia expression in wild type embryos is en- 
hanced by exogenous RA, it is not lost upon treatment 
with lOuM DEAB. 

In hoxblb'''^'^^^'''^'^^ mutants, DEAB treatment blocks 
r5 krox20 expression - similar to the effect in wild type 
embryos - while r3 size is reduced somewhat (Figure 4D). 
However, while DEAB treatment has no effect on hoxbia 
expression in wild type embryos, it blocks expression in 
hoxbib mutants (Figure 4D). In contrast, RA treatment 
has the same effect in hoxbib mutants as in wild type em- 
bryos, in that krox20 expression is lost in r3 and reduced 
in r5, while hoxbia expression is expanded (Figure 4F). 
Taken together, our results reveal that hoxbia expres- 
sion is abolished upon simultaneous removal of hoxbib 
and RA, but not when either factor is removed by itself, 
suggesting that RA and hoxbib independently activate 
hoxbia expression in zebrafish r4. 
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Figure 4 Retinoic acid and hoxbib act independently to 
activate hoxbia transcription. 19hpf wild type (A, C, E) or 
hoxbib'""' (B, D, F) embryos were treated with DMSO (control; A, B), 
lOuM DEAB (C D) or 100 nM RA (E, F) and assayed by in situ 
hybridization for expression of hoxbia in r4 (blue stain in A-F) and 
krox20 in r3/r5 (red stain in A-F). All embryos are flat mounted in 
dorsal view with anterior to the top. 



hoxbib affects nucleosome positioning at tiie promoter 
of hoxbia 

We recently found that nucleosome positioning at the 
promoter regions of zebrafish hox genes is a progressive 
process that occurs over several stages of embryogenesis 
independent of RA signaling [46]. Since Hoxbib regu- 
lates hoxbia expression, at least in part, we wanted to 
test if Hoxbib plays a role in nucleosome positioning at 
the hoxbia promoter. We took advantage of the fact that 
hoxblb'"^^^^'''''^^^^ mutant fish are viable and mapped 
nucleosomes at the hoxbia promoter of WT and hoxbib 



mutant embryos using a modified nucleosome scanning 
approach (Figure 5 A; [47]). 

We detect three well-positioned nucleosomes at ap- 
proximately 360 bp, 180 bp and 35 bp upstream of the 
transcription start site (TSS; -3, -2, and -1 nucleosomes; 
Figure 5B), as well as two less well-defined nucleosomes at 
50 bp and 200 bp downstream of the TSS (+1 and +2 
nucleosomes) at 4hpf - a stage when Hoxbib protein is 
not yet expressed and cannot yet be detected at the hoxbia 
promoter [48]. The nucleosome profile at 9hpf - when 
Hoxbib is present at the hoxbia promoter [48] - is similar 



Weicksel et al. BMC Developmental Biology 2014, 14:25 
http://www.bionnedcentral.conn/1471-213X/14/25 



Page 10 of 15 



Hox sites 



-500 -400 -300 -200 -100 0 100 200 300 




-9hpf WT 

Figure 5 Hoxbib affects nucleosome organization at the 
hoxbia promoter. A. Diagram illustrating location of primers used 
for nucleosome scanning across the hoxbia locus. B-D. Overlay of 
nucleosome profiles at the hoxblo promoter for 4hpf versus 9hpf 
wild type embryos (B), 4hpf wild type versus 9hpf hoxblb~^~ 
embryos (C) and 9hpf wild type versus 9hpf hoxblb~^~ embryos (D). 
The -3, -2, -1, +1 and +2 nucleosomes are labeled in (B). Green 
circles indicate locations where changes in nucleosome organization 
are statistically significant (p < 0.05) while arrowheads point to Site 1 
(green arrowhead) and Site 2 (red arrowhead) discussed in the text. 
Error bars represent standard error of three biological replicates. Red 
region on horizontal axis indicates approximate location of Pbx/ 
Meis/Hox binding sites. 



to that at 4hpf, but there are some notable differences 
(Figure 5B). In particular, there is reduced nucleosome 
density surrounding the TSS from approximately pos- 
ition -200 bp to position +200 bp, indicating that changes 
in nucleosome density occur at the hoxbia promoter at 
developmental time points when Hoxbib is bound at the 
promoter. 



To determine if the change in nucleosome arrangement 
between 4hpf and 9hpf is driven by Hoxbib, we compared 
the nucleosome profile of hoxbib mutant embryos at 9hpf 
to that of wild type embryos at 4hpf (Figure 5C). If 
Hoxbib is solely responsible for promoting changes in 
nucleosome arrangement, the two profiles should be 
similar. Indeed, the 9hpf hoxbib mutant and 4hpf wild 
type profiles are more similar to each other than are the 
wild type 9hpf and wild type 4hpf profiles, with the ex- 
ception of a region centered at approximately +120 bp 
downstream of the TSS (Site 2; red arrow). This finding 
suggests that Hoxbib is responsible for affecting nu- 
cleosome organization at the TSS and immediately up- 
stream, but not further downstream at Site 2. Such a 
scenario is consistent with the fact that Hoxbib is 
known to bind regulatory elements located 250-150 bp 
upstream of the hoxbia TSS. 

Accordingly, when we compare the nucleosome profile 
of 9hpf hoxbib mutants to that of 9hpf wild type em- 
bryos, we find that the region around +120 bp (Site 2) is 
relatively unaffected, again consistent with nucleosomes 
in this region being organized independently of Hoxbib 
(Figure 5D). As expected, there remains a significant 
difference upstream of the TSS (Site 1, green arrow) be- 
tween 9hpf wild type and 9hpf mutant embryos, con- 
firming that nucleosome occupancy in this region is 
likely dependent on Hoxbib binding. Together, these 
data identify two regions where nucleosome organization 
is affected at the hoxbia promoter. The first, -120 bp up- 
stream of the hoxbia TSS (Site 1), near the known 
Hoxbib binding site appears to be directly affected by 
Hoxbib binding, while the second, -120 bp downstream 
of the TSS (Site 2), appears to be independent of Hoxbib 
binding. 

Discussion 

The phenotypes of mouse Hoxal and Hoxbl germ line 
mutants, as well as the phenotypes of anti-sense medi- 
ated knock-down of zebrafish hoxbib and hoxbia have 
been characterized [20-24,26-29,43]. Differences between 
the zebrafish and mouse phenotypes suggest that these 
genes may play different roles in different species, or, 
possibly, that zebrafish phenotypes induced by anti- 
sense morpholinos (MOs) may not represent complete 
loss of function. To clarif)^ these differences we created 
targeted germ line mutations for zebrafish hoxbib and 
hoxbia with ZFN and TALEN systems, respectively. Our 
findings indicate that Hoxal and hoxbib share roles in 
hindbrain segmentation and that Hoxbl and hoxbia 
have similar roles in facial motor neuron migration. 
Comparing the phenotypes of our germ line mutants to 
those of MO loss-of-function suggests that the MO phe- 
notypes do not represent complete loss of function in all 
respects, although we cannot completely rule out the 
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possibility that differences in genetic backgrounds may 
also affect the observed differences. Our experiments 
also reveal that hoxblb and hoxbla have species-specific 
functions. In addition, we report that hoxbla transcrip- 
tion is independently regulated by both RA and Hoxblb. 
Lastly, we demonstrate that Hoxblb mediates its effect 
on hoxbla expression, at least in part, by affecting nu- 
cleosome positioning at the hoxbla promoter. 

Hoxal and hoxblb have universal as well as 
species-specific roles in hindbrain segmentation 

Comparison of the previously reported hindbrain pheno- 
types of mouse Hoxal (produced by targeted germ line 
disruption) and zebrafish hoxblb (produced by anti- 
sense MO) loss-of-function studies reveal several simi- 
larities. In particular, r3 is expanded while r4 and r5 are 
reduced in both species [20-22,24,28,29]. However, there 
are differences in the segmentation defects observed. 
In particular, the r5 segmentation defect in mouse 
Hoxal~^~ mutants is more severe, with some embryos 
losing r5 entirely [22], while hoxblbMO zebrafish em- 
bryos display a reduction, but not a loss, of r5 [29]. In 
hoxblbMO embryos r6 is also reduced, a phenotype 
not observed in Hoxal~^~ mice. The hoxblb germ line 
mutants presented here have an expanded r3 and a reduced 
r4, similar to mouse Hoxal mutants and hoxblbMO 
embryos, indicating that Hoxal and hoxblb share uni- 
versal roles in the formation of these rhombomeres. 
However, germ line disruption of hoxblb produces a 
more severe phenotype than MO injection - e.g. Mauthner 
neurons are affected in the germ line mutant, but not by 
Hoxblb MO - suggesting that the Hoxblb MO does not 
completely block Hoxblb function. Lastly, we observe a 
fully formed r5 and a reduced r6 in hoxblb germ line mu- 
tants, distinct from mouse Hoxal mutants, suggesting that 
the loss of r5 in HoxaF'^ mice is species specific and that 
hoxblb has a species-specific role in r6 of zebrafish. We 
note that a zebrafish hoxblb TILLING mutant was recently 
published [49]. While the effects on anteroposterior pat- 
terning were not characterized in detail for this mutant, 
Mauthner neuron formation was affected and it appears 
that r4 may be reduced in size - consistent with the pheno- 
type of our hoxblb ZFN alleles. 

Hoxbl and hoxbla are required for r4 specification 

Previous reports indicate that mouse Hoxbl and zebra- 
fish hoxbla share functions important for the migration 
of facial motor neurons (FMNs) from r4 during verte- 
brate hindbrain development [26,27,29] and we confirm 
this function in the hoxbla germ line mutants reported 
here. However, it has not been clear if hoxbla is abso- 
lutely required for r4 specification in zebrafish. In 
particular, r4 specification can be tested further in zebra- 
fish by assaying the formation of Mauthner neurons. In 



previous work, embryos injected with hoxblaMO were 
found to have normal Mauthner neurons (McClintock 
et al, 2002) - suggesting either that hoxbla is not abso- 
lutely required in r4, or that the MO-injections do not 
produce a complete null phenotype. Strikingly, we find 
that Mauthner neurons fail to form in hoxbla germ line 
mutant embryos, indicating that the hoxblaMO does 
not completely eliminate hoxbla activity. It is not clear 
what step in Mauthner neuron differentiation is affected 
by loss of hoxbla function, but the absence of Mauthner 
neurons was accompanied by the appearance of smaller 
neurons at the equivalent position in embryos co-injected 
with hoxblaMO and hoxblbMO [29], potentially suggest- 
ing a fate change or a failure to differentiate. Since our 
hoxbla germ line mutant embryos also lack expression of 
the r4-specific fgf3 gene, we conclude that hoxbla is re- 
quired for r4-specific gene expression and neuronal differ- 
entiation in zebrafish. 

Regulation of hoxbla transcription by RA and Hoxblb 

Initiation of Hoxbl I hoxbla transcription can be medi- 
ated by Hoxal/Hoxblb binding with Pbx and Meis/Prep 
cofactors at a Pbx/Prep/Hox responsive element (r4-regu- 
latory element) located upstream of the Hoxbl/hoxbla 
TSS [38,50-56]. In addition, expression of mouse Hoxbl 
depends on a RA response element (RARE) located 3' of 
the gene [43,57], but it is unclear if such an element is 
functional at the zebrafish hoxbla locus [29,58], suggest- 
ing that hoxbla might instead be highly dependent on 
hoxblb for its expression in zebrafish. While we find that 
the hoxbla expression domain is smaller in hoxblb mu- 
tants, the level of expression does not appear noticeably 
affected, demonstrating that hoxbla expression can be ac- 
tivated independently of Hoxblb in zebrafish. Treatment 
with DEAB revealed that the hoxbla expression observed 
in hoxblb mutants requires RA signaling. Hence, our data 
suggest that Hoxblb and RA can act independently to ac- 
tivate hoxbla expression and that they are the predomin- 
ant factors involved in this process - since simultaneous 
removal of both factors abolishes hoxbla transcription. 
However, it remains possible that some aspects of hoxbla 
expression require cooperation of RA and hoxblb - e.g. 
RA induces hoxbla more broadly in wild type than in 
hoxblb mutant embryos. Since a functional RARE has not 
been identified at the hoxbla locus, it is not clear whether 
RA acts directly, or if it acts indirectly via an intermediate 
factor. 

As discussed, Hoxblb regulates hoxbla transcription 
by forming a complex with Pbx and Prep/Meis factors at 
the r4 regulatory element upstream of the hoxbla TSS 
[48,54]. We recently demonstrated that nucleosome posi- 
tioning around the promoters of hox genes is influenced 
by DNA-binding factors during early embryogenesis, but 
that RA signaling does not appear to affect this process 
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[46]. We therefore took advantage of the fact that homo- 
zygous hoxblb'''^^^^^'"^^^^ mutant fish are viable to test if 
Hoxblb affects nucleosome organization at the hoxbla 
promoter. We find high nucleosome density across the 
hoxbla TSS at 4hpf (when Hoxblb is not yet bound to the 
promoter and hoxbla is not yet expressed), but this density 
is reduced in a region extending from +200 bp to -200 bp 
in wild type embryos at 9hpf (when Hoxblb is bound to 
the promoter and hoxbla is expressed). Strikingly, the af- 
fected region upstream of the hoxbla TSS is near the r4 
regulatory element (that is known to bind Hoxblb) and nu- 
cleosome density in this region is elevated in hoxblb mu- 
tants at 9hpf, suggesting that Hoxblb is involved in 
depleting nucleosomes in this region coincident with onset 
of hoxbla transcription. In contrast, a region downstream 
of the TSS undergoes a reduction in nucleosome density by 
9hpf in both hoxblb mutant and wild type embryos, sug- 
gesting that nucleosome changes in this region are inde- 
pendent of hoxblb and must depend on other factors. 
Taken together, our findings indicate that hoxbla expres- 
sion is regulated both by Hoxblb acting, at least in part, to 
affect nucleosome organization and by RA acting via an un- 
known pathway. 

Conclusions 

We generated several novel germ line mutants for zebra- 
fish hoxbla and hoxblb. Our analyses indicate that 
mouse Hoxbl and zebrafish hoxbla have comparable 
functions, suggesting a conserved role for these genes. In 
contrast, while mouse Hoxal and zebrafish hoxblb share 
functions in the formation of r3 and r4, they differ with 
regards to r5 and r6, where Hoxal appears to control for- 
mation of r5, but not r6, in the mouse, whereas hoxblb 
regulates formation of r6, but not r5, in zebrafish. Lastly, 
our data reveal independent regulation of hoxbla expres- 
sion by retinoic acid and Hoxblb in zebrafish. 

Methods 

Fish care 

Ekkwill (EK) embryos were collected through natural 
matings and staged using morphological criteria as de- 
fined previously [59]. 

Generation of zinc finger and Tale nucleases 

Zinc finger nucleases (ZFNs) in build 1 and 2 (Table 1) 
were constructed from a single finger modular archive 
[35]. ZFNs for build 3 were designed from a library con- 
taining optimized two -finger modules for improved DNA 
recognition [36]. Nuclease assemblies for all three builds 
were completed using previously published protocols 
[35,36]. TALENs were constructed using the Golden Gate 
TALEN assembly kit (addgene: TALEN Kit #1000000024) 
following previously published protocols [31,60]. pCS2 
plasmids containing completed ZFNs and TALENs were 



linearized and in vitro transcription was performed with 
the T7 mMachine ultra kit (Ambion: AMB1345). 50- 
100 pg of mRNA encoding ZFNs or TALENs was then 
injected into wild type embryos at the one cell stage. 

Identification of germ line mutations 

Genomic DNA (gDNA) was prepared from a pool of 50 
injected embryos at 24hpf. A fragment overlapping the 
nuclease target site (a 200 bp fragment carrying a Btgl 
restriction site from exon one of hoxbla and a 300 bp 
fragment carrying a BslI restriction site from exon one 
of hoxblb; see Additional file 6: Table S3 for primer se- 
quences) was amplified and digested to determine if the 
diagnostic restriction site had been disrupted, indicating 
an active nuclease. Embryos injected with active TALENs 
and ZFNs were raised to adulthood. These mosaic foun- 
ders (FO) were outcrossed to wild type fish and pools of 
50 embryos were genotyped as above to identify germ 
line transmission. An Fl generation was raised from each 
germ line positive founder and Fl carriers were identified 
by genotyping of gDNA isolated from fin clips [61]. Indi- 
vidual Fl carriers were then sequenced to determine if a 
frameshift had occurred. Carriers with frameshift muta- 
tions were used for phenotypic analysis. 

Genotyping 

Sequencing of Fl fish carrying frame shift mutations re- 
vealed the introduction of a Btgl restriction site in the 
hoxblb"""^'^^ and hoxblb"""^'^^ lines that was used for 
subsequent genotyping. Further, hoxblb^^^^^ contains 
an insertion that was used for subsequent genotyping by 
employing an insertion-specific primer. See Additional 
file 6: Table S3 for primer information. 

In situ probes and antibody labeling 

In situ protocols were as previous published [62]. Em- 
bryos were fixed in 4% paraformaldehyde (pfa) and 
stored in 100% methanol at -20°C. In situ probes for the 
following genes were used: hoxbla [63], krox20 [64], pax2 
[65], hoxbSa [66], fgf3 [67], and hoxd4a [68] . Visualization 
was completed using colorimetric reaction using NBT/ 
BCIP or INT/BCIP in 10% polyvinyl alcohol. Embryos 
were dissected from the yolk and flat mounted in 70% 
glycerol for imaging on bridged coverslips as described 
(Zannino and Appel, 2009). Images were captured using 
a Nikon Eclipse E600 microscope equipped with Spot 
RT Color camera (model 2.2.1). Antibody labeling with 
mouse anti-Isl (39.4D5, 1:100; Developmental Studies 
Hybridoma Bank (DSHB)) and mouse anti-3A10 (1:100; 
DSHB) was visualized using the Alexa Fluro secondary 
antibody 488 goat anti-mouse (1:200; Molecular Probes). 
Embryos were fixed in 4% AB fix (4% paraformaldehyde, 
8% sucrose, Ix PBS) for 2 hours at room temperature 
(RT) or overnight at 4C. Whole-mount fluorescent labeling 
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was performed as described [69]. Embryos were imaged 
by an upright Zeiss Imager.M2 equipped with a 20x 
water immersion objective [numerical aperture (NA) = 1.0], 
mounted on a Nano-Drive, and a Perkin Elmer Ultraview 
system. Images were imported into Adobe Photoshop and 
adjustments were limited to contrast, levels, and cropping 
and were applied to the entire image. 

Micrococcal nuclease digestions and nucleosome 
identification 

Micrococcal nuclease (MNase) digestions, isolation of 
mono-nucleosome fragments, and amplification of puri- 
fied DNA was performed on 4 and 9 hour post fertilization 
embryos as previously published [46]. Nucleosomes 
were mapped to the hoxbla promoter using a reported 
nucleosome scanning protocol [47]. Briefly, chromatin 
was isolated from embryos at 4hpf and 9hpf and 
digested with MNase. Mono-nucleosome sized frag- 
ments were purified and amplified by quantitative poly- 
merase chain reaction (qPCR) using tiled primers 
spaced -50 bp apart (See Additional file 6: Table S3 for 
primer sequences) from -450 bp upstream to -230 bp 
downstream of the hoxbla TSS. qPCR was done using 
the Qiagen QuantiFast SYBR Green PGR Kit (Qiagen: 
204054) on the ABI 7900HT Sequence detection system 
in a 384 well format and analyzed using SDS software 
v2.3. The signal from three biological replicates was av- 
eraged and values were expressed as a log2 ratio of the 
MNase sample to a control sample consisting of gDNA 
randomly fragmented by sonication. Statistical analysis 
was done using Student s ^-test in Microsoft Excel with 
significance cut-off set to p = 0.05. 

Additional files 



Additional file 6: Table S3. Primers used for ger^otypir^g and 
nucleosome scanning. 
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